A** 






AUBURN UNIVERSITY 

ENGINEERING EXPERIMENT STATION 


FRACTURE ANALYSIS OF RADIAL SCIENTIFIC INSTRUMENT MODULE 
REGISTRATION FITTINGS OF THE SPACE TELESCOPE 


FINAL REPORT: January 15, 1985 - July 15, 1986 

Contract Number: NAS8-36287 

Submitted by: Auburn University 

Engineering Experiment Station 
Auburn University, Alabama 36849 


Prepared by: 


Prepared for: 


C. W. Springfield, Jr. 

Department of Civil Engineering 

Auburn University 

Auburn University, Alabama 36849 


George C. Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 35812 


September 10, 1986 


(NASA-CR- 178952) FRACTURE ANALYSIS OF 
RADIAL SCIENTIFIC INSTRUMENT MODULE 
REGISTRATION FITTINGS CF THE SEACE TELESCOPE 
Final Report, 15 Jan. 1585 - 15 Jul. 1586 
(Auburn Univ.) 81 p CSCL 20K G3/39 



N86-32773 


Unclas 
446 o 3 





■&.c. 


FRACTURE ANALYSIS OF RADIAL SCIENTIFIC INSTRUMENT MODULE 
REGISTRATION FITTINGS OF THE SPACE TELESCOPE 


FINAL REPORT: 
Contract Number: 
Submitted by: 

January 15, 1985 - July 15, 1986 •• - i} 

NAS8-36287 

Auburn University 

Engineering Experiment Station 

Auburn University, Alabama 36849 


Prepared by: 

C. W. Springfield, Jr. 

Department of Civil Engineering 

Auburn University 

Auburn University, Alabama 36849 


Prepared for: 

George C. Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 35812 



September 10, 1986 




TABLE OF CONTENTS 


INTRODUCTION 

REGISTRATION FITTING AT POINT A 

REGISTRATION FITTING AT POINT B 

REGISTRATION FITTING AT POINT C 

REMARKS 

TABLES 

FIGURES 

REFERENCES 

APPENDIX A 


APPENDIX B 



INTRODUCTION 


The space telescope contains various scientific instrument 
(SI) modules which are mounted to the Focal Plane Structure 
(FPS) in a statically determinate manner. This is accomplished 
by using three registration fittings per SI module, one 
resisting three translations, another resisting two and the 
third resisting only one. Due to thermal insulating 
requirements these fittings are complex devices composed of 
numerous pieces. The structural integrity of these fittings is 
of great importance to the safety of the orbiter transporting 
the telescope, so in addition to the stress analyses performed 
during the design of these components, fracture susceptibility 
also needs to be considered. In this work the pieces of the 
registration fittings for the Radial SI Module containing the 
Wide Field Planetary Camera were examined to determine which 
would endanger the orbiter if they fractured and what is the 
likelihood of their fracture. The latter is stated in terms of 
maximum allowable initial flaw sizes in these pieces. 

When possible, pieces of the fittings where shown to be 
fail-safe through redundancy. Primarily this was done for the 
bolts in four bolt connections. For these it was shown that the 
loss of one bolt would not destroy the connection. For pieces 
in nonredundant configurations fatigue fracture analyses were 
performed. 

In order to determine the maximum allowable initial flaw 
sizes, fatigue fracture analyses were performed using the 
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FLAGR04 computer program written by T. Hu of Rockwell 
International [1]. This program integrates the crack growth 
rate per cycle equation of Collipriest [2] for a given geometry 
and load spectrum, beginning at a specified initial flaw size 
and ending when fracture instability occurs. Given the desired 
length of service of the parts, this program was used 
iteratively to determine the maximum allowable initial flaw 
sizes. In this work the desired length of service was taken to 
be four lifetimes, with a lifetime being defined as one 
application of the load spectrum given in Table 2. 

Some of the pieces are bolts or other threaded rods. If a 
flaw exists at the root of a thread then there is some 
interaction of the two stress raising effects; i.e., the stress 
concentration due to the thread and the stress singularity (if 
Linear Elastic Fracture Mechanics methods are used) due to the 
flaw. The extent of this interaction is not yet well defined. 
In order to obtain conservative estimates of allowable initial 
flaw size for these pieces, a fracture analysis code called 
ROD, developed by C. Meyers of MSFC, which also uses the 
Collipriest crack growth rate equation and includes the 
capability of analyzing a rod with an external circumferential 
crack, was used by treating the thread depth as part of the 
flaw size. These estimates are included in this report. 

However, inasmuch as the possible fracture of threaded parts is 
a common issue in structural analysis it was decided to 
investigate the stress concentration/stress singularity 
interaction to determine appropriate methods of stress 
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and , thus , 


intensity factor calculation for these geometries, 
to be able to make more accurate crack growth predictions, not 
only for the threaded pieces of the registration fittings being 
analyzed here, but also for future fracture analyses of such 
pai*ts. The results of this investigation are reported in 
Appendix B. 

This project was divided into four Tasks. In Task I the 
identification of required fracture analyses was accomplished. 

In Task II the FLAGR04 code, furnished by MSFC, was implemented 
on Auburn University’s IBM 3033 computer. In Task III stress 
analyses needed in addition to those used in the original 
design of the fittings were performed. These were needed to 
supply suitable input data for the fracture analyses which were 
performed as Task IV. This report, however, is not divided by 
Tasks but, instead, by fittings, these being identified by 
their location at points A, B or C, as indicated in Figure 1. 
Also shown in Figure 1 are the global coordinate directions, 

r 

V; , Vs and Vs, for the structure. Forces applied in these 
directions are identified as Ai , A 2 , A 3 at point A, Bi at point 
B, etc. Table 1 shows the loads in these directions for the 
various events in the service of the telescope, and Table 2 
shows the loading spectrum used in the fatigue fracture 
analyses of parts which experience complete load reversals. 

Some parts are loaded only during the positive or negative half 
of the load cycles. Some are subjected to a pretension. For 
these the spectrum in Table 2 was modified appropriately. The 
methods of analysis used and results obtained for each piece of 
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each fitting are reported. 

REGISTRATION FITTING AT POINT A 

The various pieces of the registration fitting at point A 
which were addressed in this work are listed in Table 3. Also 
shown are the mjaterials of which each is made, the threshold 
stress intensity factor range, AKo , of each, and the fracture 
toughness , Kc , of each, as well as the manner in which their 
fracture susceptibility was treated; i.e., a piece is listed as 
fail-safe or else its maximum allowable initial flaw size is 
specified. These are discussed as the remainder of this 
section. 

Focal Plane Structure Side 

The base (679-3973) which mounts to the FPS is shown in 
Figure 2. The most likely fracture scenario identified for the 
base was the growth of a through crack located as indicated in 
the figure. The loading on this crack was assumed to be due to 
loads applied in the Vz direction, resisted equally by the two 
shear lugs. Thus, half of the load was used as a tensile 
stress on a central through crack model as given in Figure 3, 
with w = 3.0 inches, t = 0.718 inches and vg = 0.3646 x kz 

ksi. This stress, <rg , is either the stress va,°l 0 or of 
Table 2 if kz is due to Acoustic, Lift-Off or Landing loads, 
respectively. From Table 1 it is found that the largest value 
of kz is 0.77 kips which is due to Lift-Off. Steps 8 and 16 of 
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the spectrum in Table 2 give the largest stress range, then, it 
being 0.56 ksi. With this applied loading a crack length equal 
to the plate width would not be large enough to reach the 
threshold stress intensity factor range of 7 ksi-inch 0 - 5 for 
this titanium alloy. It was concluded, then, that fracture 
would not occur in the base. 

Three flexures (679-4132) are designed to transmit force 
A 3 from the ball retainer to the cover. This causes the 
bending stress in the flexures. However, loads Ai and A 2 , which 
are primarily resisted by the three radial shear slugs, cause 
axial forces in these flexures since the applied loads, the 
shear slugs and the flexures are not coplanar. This is 
illustrated in Figure 4, and the locations of the flexures in 
their V 1 -V 2 plane is also shown. For fracture the largest 
tensile stresses are the ones of interest, and this occurs in 
the flexure identified as Fb in Figure 4. It was modeled as a 
cantilever beam with its movable end restrained against 
rotation, this end being loaded transversely and axially by 
concentrated forces. This is shown in Figure 5. The bending 
stress is <?b = 19.07 x A 3 ksi, and the tensile stress due to 
axial loading is = 9.53 x Ai ksi. For this particular 
flexure the A 2 loading did not contribute to the stress. A 
fatigue fracture analysis was performed on the model shown in 
Figure 5 considering an edge crack subjected to both tension 
and bending using FLAGR04. In the load spectrum of Table 2, 
the stresses, v A , v L0 and v L t are the sum of the bending and 
tensile stresses. The FLAGR04 program then uses the correct 
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proportions of these in tension and bending. The maximum 
allowable initial flaw size for four lifetimes was found to be, 
an = 0.022 i nches . 

The loads acting on the flexures are transferred to the 
aluminum cover (679-4135). They subject the cover alternately 
to tension and bending and then to compression and reversed 
bending, as is shown in Figure 6. The bending moment is due to 
the lateral force, F, and the force reacting it which is 

o- 

provided by the internal spacer. Their lines of action are 
assumed to be separated by a distance, e - 0.903 inches. The 
axial loading results from A 3 . Neither of the fatigue 
fracture computer programs being used has the capability to 
treat a cylinder subjected to both tension and bending, so a 
fracture model as is shown in Figure 3 was used in FLAGR04 to 
represent half of the cylinder, albeit flattened into a plate. 
In this model w “ 2.75 inches, which is half of the cover 
circumference, t = 0.25 inches, the cover wall thickness at the 
point of interest, and <rg = 1.68 x F + 0.728 x A3 ksi , with F 
= 0.857 x (Ax 2 + A2 2 ) 05 kips. It was found that for a crack 
half-length of over 1.3 inches no crack propagation will occur. 

A bolt (679-5280) holds the aluminum cover in place. The 
bolt has a pretension of 8.24 kips. When Ai,A 2 ,A 3 are positive 
this bolt is subjected to an additional tension of A 3 ^ 1.2 x F 
kips, with F defined as in the preceding paragraph. When the 
negatives of these are applied, the bolt does not carry the A 3 
load, but there is still a tensile contribution due to the 1.2 
x F load. As a result every cycle of loading produces two 
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cycles of tension in the bolt. In order to simplify the 
analysis in a conservative manner, it was assumed that the A 3 
loading produces tensile stress in the bolt on its negative 
cycle also so that every cycle would experience the same 
maximum stress. The cross sectional area of the bolt is 0.1504 
square inches, yielding a minimum stress in the bolt of 54.8 
ksi and a maximum stress of 54.8 + vg ksi, with ®g = (A 3 + 1.2 
x F)/0.1504 ksi, for each cycle of loading. Thus, the terms in 
the Maximum column of Table 2 must have the 54.8 ksi prestress 
added to them, all the terms in the Minimum column are simply 
this prestress, and in the Cycles column each number is 
multiplied by two. Using this load spectrum in the ROD program 
it was determined that fracture instability would be reached at 
four lifetimes for an initial circumferential flaw size of 
0.069 inches. However, net section yielding will occur in this 
piece before fracture instability, so analyses were performed 
to determine what initial flaw size would produce net section 
yielding at four lifetimes. This flaw size was found to be 
0.035 inches. 

The bolt which holds the aluminum cover in place mates 
with an internally threaded portion of the ball lower retainer 
(679-4130-111). This is illustrated in Figure 7. The threaded 
portion experiences the same load spectrum as the bolt except 
that the tensile area of this piece is different. In this 
analysis the minimum stress was 40.6 ksi and the maximum stress 
was given by 40.6 +(A3 + 1.2 x F)/0.2029 ksi. The ROD program 
used in the analysis of the bolt does not treat internally 
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threaded pieces, so the threaded portion of the lower retainer 
was treated as a plate of width, w = 1.61 inches, which is the 
circumference at its average diameter. Its plate thickness, t 
= 0.126 inches, is the difference between the outer radius of 
the piece (0.312 inches) and the root radius of a 7/16 inch 
bolt (0.186 inches). Instead of a central through crack, a 
through edge crack was considered, the depth of this crack 
being the thread depth plus an initial flaw depth, and the 
applied stress was assumed to be uniform. In this way it is 
felt that an approximation to an internal circumferential flaw 
was achieved. The results of this analysis predict an 
allowable initial flaw size of 0.011 inches. 

The ball upper retainer is attached to the lower retainer 
by four bolts (NAS 1351) which were checked for redundancy. 
Figure 7 shows the retainer and the location of these bolts. 

In order to demonstrate the redundancy of these bolts it was 
assumed that one of the bolts was missing and that the other 
three would carry the tension and compression required to hold 
the fitting. The missing bolt was assumed to be the one in the 
fourth quadrant of the V 1 -V 2 plane, and the Lift-Off values of 
Ai,A 2 ,As were used to compute the maximum tensile stress in a 
remaining bolt. A pretension of 2.9 kips was also applied. 

The maximum tensile stress in a bolt thus calculated was found 
to be 76.4 ksi which is less than the ultimate tensile strength 
of the bolt, Sut = 80 ksi. Consequently, it was determined 
that three bolts are capable of carrying the load, making this 
connection fail-safe, and that no fracture analysis is 
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necessary for these. 

Four bolts (NAS 1005) are used to attach the base to the 
FPS. These were checked for redundancy in the same manner as 
the ball retainer connecting bolts. The location of these 
bolts is shown in Figure 2. Loads Ai and kz were taken to be 
acting in a plane located 1.25 inches above the base/FPS 
interface for moment calculations. The bolt assumed to be 
missing is the one located in the third quadrant of the V 1 -V 2 
plane. Again using the Lift-Off values of the applied loads 
and a pretension of 3.09 kips the maximum tensile stress in a 
bolt was found to be 102 ksi which is less than the ultimate 
tensile strength of 140 ksi. Thus, this connection is 
fail-safe, and a fracture analysis of these bolts is not 
required. 

Scientific Instrument Side 

The base on the SI side of the point A fitting (679-2152) 
is shown in Figure 8. A possible fracture because of a through 
crack located as shown in the figure was investigated. The 
procedure and results are quite similar to those used and 
discovered for the base on the FPS side; that is, the loading 
was half of A 2 applied to a fracture model as shown in Figure 
3, but with w = 3.24 inches, t = 0.88 inches and <^g = 0.263 x 
A 2 ksi. As with the other base this stress is too small to 
develop a stress intensity factor range as large as the 
threshold value for any possible crack size. 

For the jackhead (679-2230) the critical location for a 
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flaw is in the thread relief grove, as indicated in Figure 9. 
This was analyzed using the ROD program assuming a 
circumferential crack with a depth equal to the groove depth 
plus an initial crack depth. Loads Ai and A 2 cause the same 
tensile stress in the jackhead during both the positive and 
negative halves of the loading cycles, while the A 3 load causes 
tensile stress during the positive half of the cycle and no 

stress during the negative half. As was done in the analysis 

I : * 

of the bolt which fastens the aluminum cover, it was ^ssumec} 
that the tension due to A 3 loading occurs in both halves of the 
loading cycle so that the number of cycles in the load spectrum 
of Table 2 may simply be multiplied by two. This results in a 
crack growth rate somewhat larger than actually exists, so a 
conservative analysis is obtained. The cross sectional area of 
the jackhead is 0.1963 square inches, and a pretension of 4.84 
kips is applied, so the Minimum stresses in the loading 
spectrum are always 24.7 ksi. The Maximum stresses in the 
spectrum are given by 24.7 + (A 3 + F)/0.1963 ksi, in which F = 

0.813 x (A 1 2 + A 2 2 ) 0 5 . The fatigue fracture analysis 
predicted a maximum allowable initial flaw depth of 0.084 
inches to reach fracture instability at four lifetimes, but as 
is the case with the aluminum cover attachment bolt, net 
section yielding will occur prior to fracture instability. In 
order to reach net section yielding not before four lifetimes 
an initial flaw depth of 0.032 inches is maximum. 

The bolts (NAS 1005) which attach the base to the SI were 
checked for redundancy in the manner used for the bolts 
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fastening the base on the FPS side of the fitting. Figure 8 
shows the locations of these bolts, and the one in the third 
quadrant of the V 1 -V 2 plane was assumed to be missing. 

Lift-Off loads were used along with a pretension of 3.07 kips. 
The Ai and A 2 loads were assumed to act in a plane 1.87 inches 
above the base/SI interface. A maximum tensile stress in a 
remaining bolt was determined to be 116 ksi which is less than 
the ultimate strength, Sut = 140 ksi, so this connection is 
also fail-safe. 

i ■ 

REGISTRATION FITTING AT POINT B 

Table 4 lists the various pieces of the registration 
fitting at point B which were considered in this work. The 
format of this table is like that of Table 3 for the fitting at 
point A. Except for the support plate on the SI side of the 
fitting, the pieces of this fitting are identical to those of 
the fitting at point C. Inasmuch as the loads are greater at 
point C, the results obtained from analyses at point C are 
taken as conservative results at point B. To see the details 
of the analyses for all the pieces of the point B fitting 
except the support plate the reader is referred to the section 
of this report titled ’’REGISTRATION FITTING AT POINT C”. The 
geometry and loading of the point B support plate (679-2228) 
are sufficiently different from the support plate at point C 
that they were analyzed independently. 

In order to identify likely fracture locations and to 
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determine the states of stress at these locations in the 


support plate at point B, a plane stress analysis of the 
support plate was performed using the SAP V finite element 
program [3]. The finite element model used is shown in Figure 
10, along with the locations at which flaws were assumed to 
exist. The cut-out region in which the flexure fits was 
modeled by reducing the Young’s modulus of the elements in that 
region by the ratio of the reduced thickness to the thickness 
of the rest of the piece. As can be seen in Figure 10, only a 
portion of the support plate was modeled, the remainder being 
treated as rigid. The most critical location found for a crack 
in this piece is indicated in the figure. The stress 
distribution at this location can be represented by that due to 
a combination of bending and axial loading, these being found 
to be given by ffb = 2.9 x Bi ksi and = 2.5 x Bi ksi. These 
only occur during half of a load cycle, so the stresses in the 
Minimum column of Table 2 were taken to be zero. An analysis 
of an edge crack was performed, and it was determined that an 
initial crack depth of 0.153 inches is acceptable. 

REGISTRATION FITTING AT POINT C 

Following the format of Tables 3 and 4, Table 5 lists the 
pieces of the point C registration fitting which were addressed 
in this project along with the material, fracture toughnesses, 
and fracture susceptibility of each. A description of the 
various analyses is given in the following paragraphs. 
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Focal Plane Structure Side 

Two possible flaw locations were investigated in the base 
(911-4236), these being illustrated in Figure 11 which shows 
two views of the base with the ball installed. At location 12 
is a through crack subjected to stresses due to the Ci loads. 
The fracture model is as shown in Figure 3 with w = 4.50 

inches, t = 0.5 inches and = 0.349 x Ci ksi. As was found 

• i 

to be the case with the other bases, even when Ci is due to 
Lift-Off this applied stress is not large enough to cause crack 
growth for any flaw size which can occur. The other flaw which 
was considered was an edge crack at location 12a. At this 
location the C 2 load causes both a uniform tension load and a 
bending load. The resulting stress is found to be <r e = 1.44 x 
C 2 ksi on the ball side of the piece, and it was assumed to 
decay linearly to zero on the back side. In the load spectrum 
of Table 2, vg is either <*a , ff Lo or , when the applied loads 
are due to either Acoustic, Lift-Off or Landing sources, 
respectively. These stresses were divided into the appropriate 
tensile and bending stresses in the FLAGR04 program in the 
analysis. It was found that an edge crack depth in excess of 
1.5 inches would be required to develop net section yielding 
which will occur before fracture instability. 

The stress relief groove on the stem of the ball 
(679-2387-110) is the most critical potential flaw location in 
this piece. A circumferential flaw was assumed to exist there, 
as is shown in Figure 12. The cyclic loading is tension due to 
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the force, F = (Ci 2 > C 2 2 ) 0 - 5 , during both the positive and 
negative halves of the loading cycles. Thus the numbers in the 
Cycles column of Table 2 were doubled, those in the Minimum 
column were the prestress of 14.4 ksi and those in the Maximum 
column were 14.4 + F/0.3068 ksi. It was determined that net 
section yielding would occur before fracture instability, at 
which time the flaw would have become 0.1485 inches deep. This 
depth is predicted to be reached at four lifetimes by a flaw of 
initial depth, ao = 0.1475 inches. 

Four bolts (NAS 1005) fasten the base to the FPS. They 
were checked for redundancy in a manner similar to those of the 
other bases. The bolt locations are shown, in Figure 11, and it 
was assumed that the Ci load acts in a plane 3.14 inches above 
the base/FPS interface. The worst condition arises when the 

bolt in the fourth quadrant of the V3-V1 plane is missing. 
Assuming this bolt to be missing and a pretension of 3.09 kips, 
the highest remaining bolt tension was determined to be 5.88 
kips, or 101 ksi. This is smaller than the ultimate tensile 
strength of the bolt, Sut = 140 ksi, so this connection is 
fail-safe. 

Scientific Instrument Side 

The ball on the FPS side fits into the support plate 
(679-2223). As was done for the support plate in the point B 
fitting, a plane stress analysis was performed using the SAP V 
finite element program. The model used is shown in Figure 13. 
As is clear from the figure only a portion of the support plate 
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was modeled, the remainder being assumed to be rigid. The 
regions which are cut-out to accept the flexures were modeled 
by reducing the Young’s modulus of the elements in these 
regions by the percentage that the material is actually 
reduced. Also shown in Figure 13 is the most critical location 
for the existence of an edge crack. At this cross section the 
stress can be represented by a contribution due to uniform 
axial stress and a contribution due to pure bending, = 1.75 
x Ci + 0.85 x C 2 ksi and ^b = 2.71 x Ci + 0.73 x Cz ksi, 
respectively. Because these only occur during half of a load 
cycle, the minimum stresses in the applied load spectrum were 
taken to be zero. Net section yielding, defined in this 
particular analysis as the development of a plastic hinge at 
this cross section, is the limiting condition here. So the 
maximum allowable initial flaw depth for an edge crack at this 
location is the depth which will grow such that the cross 
section is reduced to a size allowing net section yielding at 
four lifetimes. This initial flaw depth was found to be 0.21 
inches . 

The support plate is connected to the base (679-2211) by 
three bolts. The base attaches to the SI with four bolts. 

This is shown in Figure 14. The possible fracture due to 
though cracks emanating from a bolt hole as shown in the figure 
was considered. Loading at this bolt was assumed to be 
one-third of the applied Ci . A fracture model as is shown in 
Figure 15 was analyzed with w = 2.125 inches, t = 0.58 inches 
and the applied stress, <rg = 0.2705 x Ci . It was determined 
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that the threshold stress intensity factor range would not be 
reached for any possible initial crack size. 

Since three bolts (NAS 1005) connect the support plate to 
the base rather than four, this connection was not checked for 
redundancy, but instead a fracture analysis was performed on 
the bolt subjected to the highest loading. This bolt is 
indicated in Figure 14. In addition to the preload of 3.07 
kips, it is subjected to a fluctuating load of 0.434 x Ci kips 
during half of a loading cycle and zero during the other half. 
In view of this, the minimum stresses were taken to be the 
prestress and the maximum stresses were the sum of the 
prestress and’ the fluctuating stress. A circumferential flaw 
which is 0.027 inches deep will cause net section yielding, but 
this flaw does not grow when subjected to the stress intensity 
factor range corresponding to the applied fluctuating stresses. 
Therefore, the maximum allowable initial circumferential flaw 
depth is 0.027 inches. 

The four bolts (NAS 1005) which were used to attach the 
base to the SI were checked for redundancy. Their locations 
are indicated in Figure 14. The Ci load was assumed to act in 
a plane 2.45 inches above the base/SI interface, and the bolt 
in the first quadrant of the V 3 -V 1 plane was the one assumed to 
be missing. The largest bolt tension due to the applied loads, 
which were the Lift-Off loads, was found to be 2.1 kips. This 
load along with the preload of 3.09 kips causes a tensile 
stress in the bolt of 89.5 ksi. The ultimate tensile strength 
of the bolt is 140 ksi, so this connection is deemed fail-safe. 
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REMAHKS 


Various pieces of the registration fittings for the Radial 
SI module of the Space Telescope have been examined from a 
fracture mechanics point of view and deemed to be fail-safe or 
else have had maximum allowable flaw sizes specified for them. 
The results of these analyses are summarized in Tables 3-5 
and also in Appendix A which is comprised of tables in a form 
normally used by MSFC in summarizing fracture analysis results. 
In many instances the applied stress levels were so low that 
the threshold stress intensity factor range was never reached. 
In most of the others the allowable flaw sizes were large 
enough to be detected by visual inspection. However, for some 
parts, such as the flexures connecting the aluminum cover to 
the ball retainer in the fitting at point A, the flaw sizes 
were rather small. Eddy current, tests are capable of detecting 
flaws of this size (0.022 inches x 0.1 inches), so for those 
which have been so tested these small flaws should represent no 
danger of going undetected. 

In every instance approximations were made to err on the 
conservative side. These were pointed out in the discussions 
of the analyses for each fitting. One conservative 
approximation that was not mentioned, however, is the fact that 
retardation was not included in the crack propagation 
computations. It is probable that retardation occurs after 
Steps 8 and 16 in t. he load spectrum of Table 2, and so it is 
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expected that the predicted crack growth rates are larger than 
they are in reality resulting in smaller predicted allowable 
flaw sizes than actually may be tolerated. 



TABLES 
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TABLE 1 RADIAL SI LOADS 


Force 

Acoustic (kips) 

Lift-Off (kips) 

Landing (k 

A 1 

0.771 

2.372 

1.660 

A 2 

0.298 

0.770 

0.425 

A 3 

0.660 

2.014 

1.894 

B 1 

1.213 

3.459 

2.091 

C 1 

1.208 

3.440 

2.082 

c? 

0.958 

2.148 

0.987 
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TABLE 2 


RADIAL SI LATCHES LOAD SPECTRUM 



a. : Stresses 

A 

Calculated Using 

Acoustic Loads 



o^: Stresses 

Calculated Using 

Lift-Off Loads 



: Stresses 

Calculated Using 

Landing Loads 


Event 

Step 

Maximum 

Minimum 

Cycles 

Acoustics 

1 

1/3 x 

°A 

1/3 

x o A 

1417 


2 

2/3 x 

°A 

2/3 

X °A 

1696 


3 


°A 


d A 

487 

Ship 

4 

.39 x 

°L 

.20 

x 

155 


5 

.37 x 

°L 

.22 

x a L 

799 


6 

.35 x 

°L ' 

.24 

x a L 

13837 


7 

.33 x 

°L 

.25 

x o L 

218378 

Launch 

8 


°L0 


a L0 

4 


9 

.75 x 

°L0 

.75 

X °L0 

7 


10 

.50 x 

°L0 

.50 

X a L0 

13 


11 

.25 x 

°L0 

.25 

X °L0 

30 

Landing 

12 


°L 


a L 

4 


13 

.75 x 

a L 

.75 

x a. 

5 


14 

.50 x 

°L 

.50 

X °L 
x a L 

10 


15 

.25 x 

a L 

.25 

4 

Launch 

16 


a L0 


°L0 

4 


17 

.75 x 

°L0 

.75 

X °L0 

7 


18 • 

.50 x 

a L0 

.50 

X °L0 

13 


19 

.25 x 

a L0 

.25 
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Table 3. Point A Fitting Fracture Susceptibility 
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Table 5. Point C Fitting Fracture Susceptibility 



24 


FIGURES 


25 












wuuwuV 



Figure 5. Flexure (679-4132) and Edge Crack Model 
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Figure 6. Section Through Aluminum Cover (679-4135) 
Showing Applied Loads and Reactions 
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Figure 8. Point A Base, SI Side (679-2152) 
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Support Plate 679-2228 t = 0.31 TI 6AL-4V 0.153 1.0 TEB 4 18.7 6 ' 9 3 ' 



Circumferential Flaw These results obtained from analyses of point C fitting. 
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NOTES: *Size used in fracture mechanics model, t x w, diameter, etc. 

TC - Through Center F-S - Fail-Safe Analysis 

TEB- Through Edge Beam 
C - Circumferential Flaw 







































APPENDIX B - INVESTIGATION OF STRESS CONCENTRATION FACTOR - 
STRESS INTENSITY FACTOR INTERACTION FOR FLAWS IN FILLETED RODS 

During the coarse of this project several filleted rods and threaded 
fasteners were analyzed to establish a maximum allowable flaw size at the 
fillet or thread root for each. Insufficient work has been reported for 
stress intensity factor (SIF) solutions of these geometries to allow 
fatigue fracture analysis methods to be developed. Recently, Nord has used 
finite elements to analyze some threaded rod geometries [B-l ], and he is 
currently incorporating his results in the Del -West fatigue fracture 
program for use by MSFC. Since there is very little data with which' to 
compare his results, his results cannot be verified at this time. For the 
analyses of threaded and filleted rods in this work approximate means for 
determining the SIF were used. There were two approaches which were 
considered. One was to determine the stress concentration factor for the 
fillet and increase the applied stress in the fatigue fracture analysis by 
this factor. The fatigue fracture analysis was to be performed for a 
nonfilleted bar with an initially circular surface flaw. It was thought 
that this would lead to conservative results; i.e., prediction of crack 
growth more rapid than actually would occur. The other approach was to 
assume that in the fillet there was a circumferential surface flaw 
extending around the entire circumference. By treating the fillet depth as 
part of the flaw size a fatigue fracture analysis could be performed since 
there are solutions available for this problem. Again, it was felt that 
this would yield conservative results. The results reported for such 
pieces in this document were obtained in this latter manner. 
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Since a lack of SIF solutions for surface flaws occuring in fillets 
prevents accurate estimates of growth for these flaws, it was decided to 
investigate these flaw geometries using finite element techniques. 

Analyses of rods having circular surface flaws were analyzed to determine 
the interaction of the stress concentration factors of the fillets with the 
stress intensity factor distributions of the flaws. Figure B-l shows the 
generic filleted rod geometry which was analyzed. For all the rods 
considered the difference between the large diameter and the small diameter 
was twice the fillet radius. All dimensions were normalized by the small 
diameter. Thus, there were two dimensionless parameters which 
characterized the rod geometry; r/d, and a/d. The lengths of the rods were 
long enough to not affect the results. Typical elastic constant for steel 
(E = 30,000,000 psi, v = 0.3) were used. 

Stress intensity factor distributions around the crack fronts for the 
various geometries were determined from three-dimensional finite element 
analyses. By assuming two surface flaws in the bars located 180° apart the 
various symmetries allowed a model of one eighth of a rod to be used. The 
sizes of the flaws considered were small enough that no interaction between 
the flaws should have occured, so it is considered that the results of only 
one flaw in a rod would be essentially the same. A typical finite element 
mesh used is shown in Figure B-2. The region surrounding the crack front 
was modeled with collapsed sided, 20 noded, isoparametric elements with the 
mid-edge nodes of the radial edges moved to the quarter point to produce 
the proper strain singularity at the crack front. This modeling technique 
was suggested by Barsoum [B-2 ] . Inasmuch as this is a displacement based 
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finite element method, the SIF must be determined by post-processing the 
displacements in some way. The SIF were determined from these results in 
two ways; by using the nodal displacements in the analytical expressions 
for displacement, and by using an energy release rate method developed by 
deLorenzi [B-3, B-4]. Each of these methods is discussed in the following 
paragraphs. 

The displacement functions in the finite elements are qssumed in terms 
of their values at their nodes, which are the displacement values 
determined in the finite element solution. When the modeling technique of 

Barsoum is used the displacement function along a radial line originating 

' ' i ; i : l » 

at the crack front is given by 

u = ui +.[4U 2 - u 3 - 3 u i]V-f + [2 u 3 + 2«! - 4 u 2 ] -f 

In this, ui is the displacement at the crack front node, U 2 is the 
displacement at the quarter point node, and U 3 is the displacement at the 
corner node away from the crack front. On the other hand the expression 
for displacement in the vicinity of a crack is known to be of the form, 



The parameter, C, depends on the proximity to the free surface along the 
crack front. Equating the coefficients of the r terms in the above 
equations yields: 

Kj = JH [4u2 - U3 - 3ui] 

C-l ? L 

Thus, one way to determine the SIF is to utilize the displacements 
determined by the finite element solution in this equation. This method is 
sometimes referred to as the "displacement substitution method" and will be 
in this report. 
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Since the stress intensity factor is related to the amount of energy 
released when a crack grows an incremental amount, calculating the energy 
release rate for a crack provides a reliable means of determining SIF. 
Several schemes for calculating the energy release rate for flaws have been 
reported over the years. Most of these, in one way or another, calculate 
the difference in strain energy of two bodies identical in every way except 
that the crack size of one is slightly larger than that of the other. 
Dividing this difference in strain energy by the difference in the crack 
sizes yields an approximation of the energy release rate. In the limit, as 
the difference in crack sizes approaches zero, this approximation 
approaches the exact energy release rate. Thus, the accuracy of the 
calculation depends on the increment in crack size used. Recently, 
utilizing the fact that the energy release rate .is given identically by the 
J - integral, deLorenzi [B-3, B-4 ] has developed a method for calculating 
energy release rate rather simply. What is more, the increment in crack 
size divides out of his equations so that the results do hot depend on the 
magnitude of this quantity at all. This method will be called the 
"deLorenzi method". 

Before the filleted rods were analyzed, both the displacement 
substitution method and the deLorenzi method were used to analyze 
semicircular surface flaws of a/d = 1/12 in nonfilleted rods, the results 
being compared with each other and with results of other investigators 
[B-5, B-6 ] to assess the accuracy of the two methods. This comparison is 
displayed in Figure B-3. In the figure the SIF distributions around the 
crack front are displayed in non-dimensional form as geometry factors, F. 
The angle, 0, is measured from the tangent to the rod surface where it 
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intersects the origin of the semicircular flaw radius, as shown in Figure 
B-l. Two finite element meshes are compared, one using 102 elements and 
the other 204 elements. The various curves fall in a band which has a 
Width of approximately 5% of the mean value. The mean value corresponds 
well with the analytical solution of a semicircular surface flaw in a 

semi -inf ini te body [b- 6 ] and also with the results of [B-5 ]. The lower 

\ ’ 

edge of the band is the curve corresponding to the 102 element mesh, 
deLorenzi method results, while the upper edge is the curve from [B-5 ]. If 
one takes the curves from [B-5, B-6] as essentially correct and if one 
bears in mind the fact that the results in [B-5] were obtained using the 
deLorenzi method with a 340 element mesh, then it is seen that as the mesh 
is refined when using the deLorenzi method the solution converges toward 
the correct one. This is not the case using the displacement substitution 
method. Moreover, the displacement substitution results oscillate in the 
region near the rod surface. Consequently, it is concluded that the 
deLorenzi method is the more reliable one. Also shown in Figure B-3 are 
curves for a rod with a fillet of r/d = 1/4 as determined using the 
deLorenzi method. It is clear that the 104 element and 208 element curves 
have the same relative locations as for the nonfilleted rod. In view of 
this it was decided that for the various filleted rod analyses the results 
obtained using the 102 element grid and the deLorenzi method could be 
considered to be approximately 5 % low but otherwise correct. Although the 
results obtained using the displacement substitution method are included in 
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this appendix, because of their oscillatory nature it is thought that the 
results obtained with the deLorenzi method are more reliable. Therefore, 
it is recommended that for analysis purposes one use the curves reported as 
determined by the deLorenzi method and multiply them by 1.05. 

Four different flaw sizes, a/d = 1/6, 1/8, 1/12, 1/20, were 
investigated along with three different fillet sizes, r/d = 1/4, 1/7, 1/10, 
as well as the nonfilleted rod condition, r/d = 0. It is generally thought 
that the stress intensity factors for surface flaws of size a/d < 1/10 can 
well approximate solutions for flaws in semi-infinite bodies, and the 
results given in Figure B-3 from references [8-5 ] and [B-6 ] support this. 
Figures B-4 and B-5 show the results determined in this study for 
nonfilleted rods, and they support this conclusion also except for the 
curves corresponding to a/d = 1/20. The curves for the other flaw sizes 
appear to be converging to some lower bound curve, possibly that of a/d = 
1/12, but the a/d = 1/20 curve is far below this. It is felt that this 
curve is incorrect, but the reason has not been determined at this time. 

Figures B-6 through B-9 are the results determined for the filleted rods 
using the deLorenzi method. The displacement substitution method 
counterparts are given in Figures B-10 through B-13. In each of these it 
is seen that the SIF values where the flaws intersect the fillet surfaces 
are much higher than for the nonfilleted rods, but not by an amount 
equivalent to the stress concentration factor for each particular fillet. 
Furthermore, it is seen that at the deepest penetrations of the flaws into 
the rods the SIF values are much lower, and, in fact, for fillet sizes 
approximately equal to or less than the flaw sizes the SIF's are actually 
smaller than those in nonfilleted rods. At first this seems a surprising 
result, but after some consideration it becomes clear why this is so and 
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why the values near the surfaces are lower than simply the nonfilleted rod 
values multiplied by the appropriate stress concentration factors. 

An explanation of this behavior is best given after first recalling 
that the stress intensity factors of a flaw in a body loaded remotely by 
some stress distribution, o 0 , are exactly the same as those obtained if one 
instead loads the flaw surface with tractions equal to the negative of 
those caused in an unflawed body by o 0 . Specifically, one could obtain 
the S IF values of the unfilleted rods by applying a uniform pressure 
distribution on the flaw surfaces just as well as by applying a remote 
uniform tension, since in unflawed, nonfilleted rods the remote loading is 
transmitted undisturbed throughout the rods. Thus, one can consider that 
the SIF results given in Figure B-3, obtained in any way by anyone, are for 
flaws loaded on their surfaces by uniform pressures equal to the stress in 
Kj = - a oN /^. This is illustrated in Figure B-14a. Similarly, if one 

wished to determine SIF values for flaws which occur at the roots of 
fillets in rods loaded in remote tension, one could do so by applying the 
negative of the stress distribution caused by that tension at the flaw 
location. In this case, however, the distribution would not be uniform, as 
is shown in Figure B-14b. It would have a maximum value at the fillet 
surface (equal to the stress concentration factor times the references 
stress) and would decrease at locations closer to the center of the rod. 
Comparing Figures B-14a,b it is clear that to simply multiply the SIF 
distributions for nonfilleted rods by the stress concentration factor would 
be to apply too much load on the surface and obtain considerable over- 
estimates. On the other hand if one uses the value of stress at the 
deepest penetration of the flaw one is using too little stress; however, in 
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some instances the value of the SIF at that location is predicted fairly 
well anyway. 

This is also a surprising result at first glance, but is understandable 
if one recalls (see, for example, Broek [B-7 ], pages 78-80) that the stress 
near the location on the crack front where the SIF is being determined has 
more effect than stress at some distance from that location. So although 
the stress increases at locations closer to the rod surface, the effect on 
the SIF at the deepest penetration decreases, and the result is very 
similar to that of a uniform stress of the smallest magnitude acting on the 
flaw. For fillet sizes approximately equal to or less than the flaw sizes 
the value of stress at this location is less than the remotely applied 
tension, so the SIF is less than would occur in the nonfilleted rod. 
Likewise, the highest stress values have the most effect on the SIF near 
the surface and the lower stress values have less effect. The curves for 
filleted rods in Figures B-6 through B-9 and B-10 through B-13 demonstrate 
these facts very clearly, so they are not surprising after all. 

Figures B-15 through B-17 and B-18 through B-20 illustrate the same 
results as in Figures B-6 through B-9 and B-10 through B-13 but grouped by 
fillet size rather than flaw size. Grouping them in this way shows clearly 
that the geometry factors decrease with increasing flaw size. This is due 
to the fact that the smaller the flaw size the more the concentrated stress 
at the surface affects its SIF distribution. Again, note that the a/d = 
1/20 results are questionable, but the others are considered to be 
reliable. 

The stress intensity factor distribution results presented in this 
appendix are useful for calculating SIF's for a few combinations of 
semicircular surface flaws occuring in filleted rods and for those 
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geometries which can be obtained by interpolation of these. More 
significantly, however, they illustrate the importance of the stress 
distribution in the region where the flaw is located. In order to 
determine SIF values for flaws in regions of stress concentration, then, it 
is best to determine the stress distribution in that region without a flaw 
and then use that distribution as the applied loading on the flaw surface. 
If a Green's function can be determined for the flaw geometry of interest, 
then any distribution of loading on the flaw surface can be handled with 
relatively little expense. Such functions have been determined for only a 
limited number of surface flaw conditions to date, however, so this method 
of SIF determination for surface flaws is not currently available for 
routine utilization. The Schwartz alternating method was demonstrated to 
be useful in determining SIF's for elliptical surface flaws by Shah and 
Kobayashi [ B-8 ] , and it has since been refined for such flaws by 
Vijayakumar and Atluri [B-9 ] and Nishioka and Atluri [B-IO], In view of 
the results shown in this report, it should be particularly powerful for 
elliptical surface flaws pccuring in regions of stress concentration since 
it utilizes an analytical solution for a buried elliptical flaw subject to 
arbitrary loading on its surfaces along with a finite element solution of 
the actual geometry without a flaw. In order to improve fatigue crack 
growth analysis capabilities a relatively inexpensive method for 
determining SIF values for surface flaws of arbitrary shape located in 
regions of arbitrary stress distribution needs to be devised, and it is 
recommended that this be pursued in future research efforts. 
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Figure B-3. Analysis Comparisons 
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Figure B-5. Geometry Factor for Nonfilleted Rod by Displacement 
Substitution Method 






Figure B-6. Geometry Factor for a/d = 1/6 Flaw by de Lorenzi Method 
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Figure B-8. Geometry Factor for a/d = 1/12 Flaw by de Lorenzi Method 
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Figure B-9. Geometry Factor for a/d = 1/20 Flaw by de Lorenzi Method 




DJI 



0 (degrees) 


Figure B- 10. Geometry Factor for a/d = 1/6 Flaw by Displacement Sub- 
stitution Method 
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Figure B-12. Geometry Factor for a/d = 1/12 Flaw by Displacement Sub 
stitution Method 
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Figure B-13. Geometry Factor for a/d = 1/20 Flaw by Displacement Sub- 
stitution Method 
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Figure B- 15. Geometry Factor for r/d = 1/4 Filleted Rod by de Lorenzi 
Method 
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Figure B- 16. Geometry Factor for r/d = 1/7 Filleted Rod by de Lorenzi 
Method 
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Figure B-18. Geometry Factor for r/d = 1/4 Filleted Rod by Displacement 
Substitution Method 
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Figure B-19. Geometry Factor for r/d = 1/7 Filleted Rod by Displacement 
Substitution Method 
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Figure B-20. Geometry Factor for r/d = 1/10 Filleted Rod by Displacement 
Substitution Method 
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